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Abstract 

We study the Bose-Einstein condensate (BEC) in a 3-D toroidal Mexican hat 
trap. By changing the parameters of the potential, or the number of bosons, it 
is possible to modify strongly the density profile of the BEC. We consider the 
ground state properties for positive and negative scattering length and cal- 
culate the spectrum elementary excitations. We also discuss the macroscopic 
phase coherence and superfluidity of the BEC by analyzing vortex states and 

their stability. 
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I. INTRODUCTION 



Recent spectacular experiments with alkali vapors 87 Rb, 23 Na and 7 Li confined in mag- 
netic traps and cooled down to a temperature of the order of 100 nK [1,2,3] have renewed 
the interest in the Bose-Einstein condensation. Theoretical studies of the Bose-Einstein 
condensate (BEC) in harmonic traps have been performed for the ground state [4,5,6,7], 
collective low-energy surface excitations [8,9] and vortex states [6]. The presence of vortex 
states is a signature of the macroscopic phase coherence of the system (the existence of a 
macroscopic quantum phase has been recently demonstrated [10]). Moreover, vortices are 
important to characterize the superfluid properties of Bose systems [11]. It has been found 
that the BEC in monotonically increasing potentials can not support stable vortices in the 
absence of an externally imposed rotation [12]. Instead, stable vortices of Bose condensates 
can be obtained in 1-D [13] and quasi-2D [14] toroidal traps: such Bose condensates are 
superfluid [11]. 

In this paper we study a 3-D toroidal trap given by a quartic Mexican hat potential along 
the cylindrical radius and a harmonic potential along the z axis. The resulting trapping 
potential is very flexible and it is possible to modify considerably the density profile of the 
BEC by changing the parameters of the potential or the number of bosons. We analyze 
the ground state properties and the vortex stability of the condensate for both positive and 
negative scattering length and calculate also the spectrum of the Bogoliubov elementary 
excitations. In particular, we consider 87 Rb and 7 Li atoms. 

The Gross-Pitaevskii energy functional [15] of the Bose-Einstein Condensate (BEC) 
reads: 

I = I dh l^ |w(r)|2 + WI*WI 2 + IT^MI 4 ' (1) 

where \l/(r) is the wave function of the condensate normalized to unity, Vo(r) is the external 
potential of the trap, and the interatomic potential is represented by a local pseudopotential 
so that g = Aixh 2 a s /m is the scattering amplitude (a s is the s-wave scattering length). N is 
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the number of bosons of the condensate and m is the atomic mass. The scattering length a s 
is supposed to be positive for 87 Rb and 23 Na, but negative for 7 Li. It means that for S7 Rb 
and 23 Na the interatomic interaction is repulsive while for 7 Li the atom-atom interaction 
is effectively attractive. The extremum condition for the energy functional gives the Gross- 
Pitaevskii (GP) equation [15] 

" - |V + V (r) + gN\^(r)\ 2 } ¥(r) = p*(r) , (2) 
L 2m J 

where p is the chemical potential. This equation has the form of a nonlinear stationary 

Schrodinger equation. 

We study the BEC in an external Mexican hat potential with cylindrical symmetry, 
which is given by 

Voir) = ^- P ir + ^z\ (3) 

where p = ^/x 2 + y 2 and z are the cylindrical coordinates. This potential is harmonic along 
the z axis and quartic along the cylindrical radius p. V (r) is minimum along the circle of 
radius p = p at z = and Vo(r) has a local maximum at the origin in the (x,y) plane. 
Small oscillations in the (x,y) plane around po have a frequency u± = po(2\/m) 1 ^ 2 . 

First, let us consider the Thomas-Fermi (TF) approximation: i.e. neglect the kinetic 
energy. It is easy to show that the kinetic energy is negligible if N » (h 2 /2m)(Xpl + 
muj 2 /2)/pl, where p = {2/ii 2 )(\/A) 1/A (muj 2 /2) 1/A g 1/2 is the bare chemical potential. This 
condition is satisfied for A 2 p(j >> lQh 2 (\pl + muj 2 /2)/ (2m). In the TF approximation we 
have 
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V2, 



* (r) = \-g~N^ ~ V ° {r)) \ 9(/i " V ° {r)) ' (4) 
where Q(x) is the step function. For our system we obtain that: a) the wave function has 
its maximum value at p = p and z = 0; b) for p < App/4 the wave function has a toroidal 
shape; c) for p > Apg/4 the wave function has a local minimum at p = z = 0; d) the chemical 
potential scales as p ~ poiV 1 / 2 . It is important to note that the TF approximation neglects 
tunneling effects: to include these processes, it is necessary to analyze the full GP problem. 



II. GROUND STATE PROPERTIES AND ELEMENTARY EXCITATIONS 

We perform the numerical minimization of the GP functional by using the steepest 
descent method [16]. It consists of projecting onto the minimum of the functional an initial 
trial state by propagating it in imaginary time. In practice one chooses a time step At and 
iterates the equation 

^(r,r + At) = *(r,r) - AtHV(t, r) , (5) 

by normalizing ^ to 1 at each iteration. 

We discretize the space with a grid of points taking advantage of the cylindrical sym- 
metry of the problem. At each time step the matrix elements entering the Hamiltonian are 
evaluated by means of finite-difference approximants. We use grids up to 200 x 200 points 
verifying that the results do not depend on the discretization parameters. The number of 
iterations in imaginary time depends on the degree of convergence required and the good- 
ness of the initial trial wave function. We found that strict convergence criteria have to be 
required on the wave function in order to obtain accurate estimates of the wavefunction. 

In our calculations we use the ^-harmonic oscillator units. We write p in units a z = 
(hKmuOz)) 1 / 2 = 1 /an, A in units {huo z )a~^ = 0.477 (5.92) peV//im 4 and the energy in 
units hu z = 0.477 (5.92) peV for 87 Rb ( 7 Li). Moreover, we use the following values for the 
scattering length: a s = 50 (-13) A for 87 Rb { 7 Li) [1,3]. 

We have to distinguish two possibilities: positive or negative scattering length. In the 
case of positive scattering length we can control the density profile of the BEC by modifying 
the parameters of the potential and also the number of particles. In Figure 1 we show the 
ground state density profile of the 87 Rb condensate for several numbers of atoms. For small 
number of particles the condensate is essentially confined along the minimum of Vb(r), there 
is a very small probability of finding particles in the center of the trap so that the system 
is effectively multiply connected. As N increases the center of the trap starts to fill up 
and the system becomes simply connected. The value of iV for which there is a crossover 
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between the two regimes increases with the value of A and of po and, within Thomas- Fermi 
approximation, scales like A 3//2 Pq. In Table 1 we show the energy per particle, the chemical 
potential and the average transverse and vertical size for the trapping potential characterized 
by the parameters po = 2 and A = 4 in the z-harmonic oscillator units. As expected, the 
energy per particle and the chemical potential grow by increasing the number of particles 
but they do not scale as N 1 ^ 2 because, with this trapping potential, the TF approximation 
is valid for iV >> 10 4 . It is instead interesting to observe that V< x 2 > = V< V 2 > grows 
less than V < z 2 > due to the presence of a steep (quartic) potential along the transverse 



direction p = \Jx 2 + y 2 and a softer (quadratic) barrier along the vertical direction z. 

In the case of negative scattering length, it is well known that for the BEC in harmonic 
potential there is a critical number of bosons N c , beyond which there is the collapse of the 
wave function [7]. We obtain the same qualitative behavior for the 7 Li condensate in our 
Mexican hat potential. However, in cylindrical symmetry, the collapse occurs along the 
line which characterizes the minima of the external potential, i.e. at p = p and z — 0. 
The numerical results are shown in Figure 2. We notice that, for a fixed p , the critical 
number of bosons N c is only weakly dependent on the height of the barrier of the Mexican 
potential. These results suggest that we can not use toroidal traps to significantly enhance 
the metastability of the BEC with negative scattering length. 

To calculate the energy and wavefunction of the elementary excitations, one must solve 
the so-called Bogoliubov-de Gennes (BdG) equations [18,19]. The BdG equations can be 
obtained from the linearized time-dependent GP equation. Namely, one can look for zero 
angular momentum solutions of the form 

(r, t) = e-i"* [p{p, z) + u(p, z)e~ iult + v*(p, z)e iu)t ] , (6) 

corresponding to small oscillations of the wavefunction around the ground state solution 
ip. By keeping terms linear in the complex functions u and v, one finds the following BdG 
equations 

-^W 2+ W P + ^) + v ^ z) -» + 2aNmp > z ^\<* z)+ 



+ gN\ip{p, z)\ 2 v(p, z) = hu u(p, z) , 
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[ - L + w P + U + v ^ z) -» + 2gN ^ ^ z)+ 



+ gN\il>(p, z)\ 2 u(p, z) = -hu v(p, z) . 



(8) 



The BdG equations allow one to calculate the eigenfrequencies uo and hence the energies 
Tiuj of the elementary excitations. This procedure is equivalent to the diagonalization of the 
N-body Hamiltonian of the system in the Bogoliubov approximation [17]. The excitations 
can be classified according to their parity with respect to the symmetry z — > —z. 

We have solved the two BdG eigenvalue equations by finite-difference discretization with 
a lattice of 40 x 40 points in the (p, z) plane. In this way, the eigenvalue problem reduces 
to the diagonalization of a 3200 x 3200 real matrix. We have tested our program in simple 
models by comparing numerical results with the analytical solution and verified that a 40 x 40 
mesh already gives reliable results for the lowest part of the spectrum. In Table 3 we show 
the lowest elementary excitations of the Bogoliubov spectrum for the ground state of the 
system. One observes the presence of an odd collective excitation at energy quite close 
to hu — 1 (in units hu z ). This mode is related to the oscillations of the center of mass 
of the condensate which, due to the harmonic confinement along the z-axis, is an exact 
eigenmode of the problem characterized by the frequency cu z , independently of the strength 
of the interaction. For large N the lowest elementary excitations saturate, suggesting that 
the Thomas-Fermi asymptotic limit is reached. 

In the case of negative scattering length we verified that, quite close to the critical number 
of bosons N c , an even mode softens driving the transition towards a collapsed state. 



Let us consider states having a vortex line along the z axis and all bosons flowing around 
it with quantized circulation. The observation of these states would be a signature of macro- 
scopic phase coherence of trapped BEC. The axially symmetric condensate wave function 



III. VORTICES AND THEIR METASTABILITY 
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can be written as 



**(r) = Mp> z ) e 



(9) 



where 6 is the angle around the z axis and k is the integer quantum number of circulation. 
The resulting GP functional (1), representing the energy per particle, can be written in 
terms of ipkij) by taking advantage of the cylindrical symmetry of the problem: 



Due to the presence of the centrifugal term, the solution of this equation for k ^ has to 
vanish on the z axis providing a signature of the vortex state. 

Vortex states are important to characterize the macroscopic quantum phase coherence 
and also superfluid properties of Bose systems [11]. It is easy to calculate the critical 
frequency fl c at which a vortex can be produced. One has to compare the energy of a vortex 
state in a frame rotating with angular frequency Q, that is E — QL Z , with the energy of the 
ground state with no vortices. Since the angular momentum per particle is hk, the critical 
frequency is given by hfl c = (E k /N — E /N)/k, where E k /N is the energy per particle of the 
vortex with quantum number k. In Table 2 we show some results for vortices of 87 Rb. The 
critical frequency turns out to increase slightly with the number of atoms. This corresponds 
to a moderate lowering of the momentum of inertia per unit mass of the condensate when 
N grows. 

For 7 Li we calculate the critical number iV c of bosons for which there is the collapse of 
the vortex wave function. We find that N c has a rather weak dependence on the quantum 
number of circulation k. Note that, in the case of a harmonic external potential, there is 
an enhancement of N c by increasing k because in that case rotation strongly reduces the 
density in the neighborhood of the origin, where the external potential has its minimum [6] . 

Once a vortex has been produced, the BEC is superfluid if the circulating flow persists, 
in a metastable state, in the absence of an externally imposed rotation [11]. As discussed 





(10) 
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previously, vortex solutions centered in harmonic traps have been found [6], but such states 
turn out to be unstable to single particle excitations out of the condensate. To study the 
metastability of the vortex we first analyze the following Hartree-Fock equation [12] 

■ " !^ V2 + V ° {r) + 2 ^l<Mr)| 2 ]0(r) = e0(r) , (11) 

which describes, in the weak coupling limit, one particle transferred from the vortex state 
^fc(r) to an orthogonal single-particle state 0(r). Quasiparticle motion is governed by an 
effective Hartree potential f e //(r) = Vo(r) + 2g\ip k (r) | 2 , which combines the effects of the 
trap with a mean repulsion by the condensate. Figure 3 shows v e ff(p, z) for N = 5000 and 
50000. The repulsion induced by the underlying condensate is quite evident near p = p . 
Let pk be the chemical potential of the vortex state characterized by a circulation quantum 
number k, then the vortex is metastable if e > p k and unstable if e < p k [12]. As shown 
in Table 2, for our 3-D system all the studied vortices are metastable and so the BEC can 
support persistent currents, thus it is superfluid. Contrary to what may be inferred by 
means of semiclassical arguments [12] the wave function describing the excitation 0(r) is not 
localized near the symmetry axis even for rather large numbers of atoms. A bound state 
at p = z = should pay a large kinetic energy cost due to the strong localization of the 
particle induced by the effective potential. Instead, it is more convenient to place the excited 
particle on top of the Bose condensate i.e. at p = p and z ^ as shown in Figure 4. 

It is well known that the Hartre-Fock approximation describes only single particle exci- 
tations [17]. To have the complete spectrum, including collective excitations, one must solve 
the BdG equations [18,19]. One must look for solutions of the form 

# fe (r, t) = e~^ kt [MP, z)e ike + u{p, z)e^ k+q ^ e^* + v*(p, z)e l{k ~ q)e e lu)t ] . (12) 

Here q represents the quantum number of circulation of the elementary excitation. We 
have solved the two BdG eigenvalue equations by finite-difference discretization using the 
same method described in Section II. We have checked that a 40 x 40 mesh gives the 
correct excitation energies within Hartree-Fock approximation. Therefore, for the purpose of 



determining the stability of the vortex state, this rather coarse mesh is sufficiently accurate. 
The results are shown in Table 4: The lowest Bogoliubov excitation is positive and always 
lower than the lowest Hartree-Fock one. Moreover by increasing the number of particles their 
difference increases as expected for collective excitations. We have also verified that vortex 
states become unstable by strongly reducing either density (down to about one hundred 
bosons in our model trap) or scattering length. 

Therefore, the behavior of the 3-D trap we have analyzed closely resembles the simplified 
1-D model studied in Ref. [13] which represents the limit of deep trapping potential. Also in 
that case Bogoliubov approximation has been used to evaluate the spectrum of elementary 
excitations showing that vortices are stabilized by strong repulsive interparticle interactions 
(or equivalently by high density). The 1-D model, however, should be taken with caution be- 
cause other branches of low energy collective excitations are present in such low- dimensional 
systems [20]. 

IV. CONCLUSIONS 

We have studied the Bose-Einstein condensate in a 3-D toroidal trap given by a quartic 
Mexican hat potential along the cylindrical radius and a harmonic potential along the z axis. 
We have shown that it is possible to modify strongly the density profile of the condensate 
by changing the parameters of the potential or the number of bosons. The properties of 
the condensate and its elementary excitations have been analyzed for both positive and 
negative scattering length by considering 87 Rb and 7 Li atoms. For 7 Li, which has negative 
scattering length, we have calculated the critical number of atoms for which there is the 
collapse of the wave function. The results have shown that a toroidal trap does not enhance 
the metastability of the ground state of the condensate. On the other hand, in the case 
of a harmonic external potential, we have recently shown [21,22] that, when a realistic non 
local (finite range) effective interaction is taken into account, a new stable branch of Bosc 
condensate appears for 7 Li at higher density. Presumably a similar state can be found also 
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in presence of a toroidal external trap for a sufficiently large number of particles when non 
locality effects are included. 

A superfluid is characterized by the presence of persistent currents in the absence of an 
externally imposed rotation. In order to investigate this peculiar sign of the macroscopic 
phase coherence of the condensate, we have also studied vortex states. Our results suggest 
that vortices can support persistent currents in 3-D toroidal traps with fairly large numbers 
of atoms. This feature essentially depends on the toroidal geometry of the trap and should 
be independent on other details of the confining potential. 
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FIGURE CAPTIONS 



Figure 1. Particle probability density in the ground state of Rb atoms as a function 
of the cylindrical radius at z = (symmetry plane). The curves correspond to different 
numbers of atoms: from 5000 to 50000. Parameters of the external potential: po = 2 and 
A = 4. Lengths are in units of a z — 1 pm and A in units of {Tiuj z )a~ A = 0.477 peV/pm 4 . 

Figure 2. Critical number N c of 7 Li atoms versus the potential barrier at the origin: 
V (0) = Apq/4. Open squares: p = 2; full squares: p = 3; open circles: p = 4. Energy is 
in units of Hcu z = 5.92 peV {uj z = 9.03 kHz) and length in units of a z — 1 pm. 

Figure 3. Effective potential v e ff(p,z) appearing in the eigenvalue equation for the 
single particle excitation Eq. (8). Two sections at z = and z = 3 (z = 6) are shown for 
iV = 5000 (N = 50000) atoms in panel a (b). z = corresponds to the symmetry plane. 
The dotted line represents the external potential. The chemical potential of the vortex state 
is marked by a short dashed line, the excitation energy by a long dashed line. Parameters 
of the external potential: po = 2 and A = 4. Units as in Fig. 1 with hu z = 0.477 peV 
(uj z = 0.729 kHz). 

Figure 4. Particle probability density of the k — 1 vortex state (solid line) and square 
of the excitation wave function (dashed line) at the radial distance p = po = 2, i.e. where 
both the wave functions peak. Curves are for N = 5000 (N = 50000) atoms in panel a (b). 
Units and parameters as in Fig. 3. 
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N E/N fi 


V< P 2 > 


V< z 2 > 


5000 5.85 7.71 


1.96 


1.41 


10000 7.45 10.26 


1.97 


1.70 


20000 9.84 14.00 


1.97 


2.05 


30000 11.73 16.94 


1.97 


2.29 


40000 13.36 19.47 


1.98 


2.48 


50000 14.81 21.74 


1.99 


2.63 



Table 1. Ground state of 8 Rb atoms in the toroidal trap with p = 2 and A = 4. 
Chemical potential and energy are in units of %uj z = 0.477 peV (uo z = 0.729 kHz). Lengths 
are in units of a z — 1 jum. 



n ejn m e nn c 

5000 6.00 7.87 9.56 0.15 

10000 7.61 10.44 12.46 0.16 

20000 10.02 14.22 16.60 0.18 

30000 11.93 17.20 19.80 0.20 

40000 13.57 19.75 22.54 0.21 

50000 15.04 22.04 25.09 0.23 

Table 2. Vortex states and excitation energies of 87 Rb atoms with k — 1 in the toroidal 
trap with po = 2 and A = 4 within Hartree-Fock approximation. Units as in Tab. 1. 
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N hui T1U2 hu> 3 hu 4 



1 1.00 

5000 1.00 
10000 1.00 
20000 1.00 
30000 1.00 
40000 1.00 
50000 1.00 



1.98 2.97 3.96 
1.70 2.43 3.19 
1.68 2.37 3.08 
1.66 2.32 3.00 
1.66 2.30 2.96 
1.66 2.30 2.95 
1.66 2.30 2.95 



Table 3. Lowest elementary excitations of the Bogoliubov spectrum for the ground state 
of 87 Rb atoms in the toroidal trap with p = 2 and A = 4. Units as in Tab. I. 



N huj e — ji\ 



5000 1.22 1.69 

10000 1.48 2.02 

20000 1.73 2.38 

30000 1.88 2.60 

40000 1.99 2.79 

50000 2.08 3.05 



Table 4. Bogoliubov vs Hartree-Fock lowest elementary excitation for a vortex state of 
87 Rb atoms with k — 1 in the toroidal trap with po = 2 and A = 4. Units as in Tab. I. 
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N=10000 
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